Collective modes of spin, density, phase and amplitude in exotic superconductors 
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The equations of motion of pair-like excitations in the superconducting state are studied for various 
types of pairing using the random phase approximation. The collective modes are computed of a 
layered electron gas described by a t — tight-binding band, where the electrons experience besides 
the long-range Coulomb repulsion an on-site Hubbard U repulsion and a nearest-neighbour attractive 
interaction. From numerical calculations we see, that the collective mode spectrum now becomes 
particularly rich. Several branches can occur below the continuum of quasi-particle excitations, 
corresponding to order-parameter fluctuations of various symmetries of pairing, and collective spin- 
density fluctuations. From the collective mode softening near the nesting vectors it is concluded, 
that in the d-wave paired state an instability occurs toward the formation of a spin-density wave. 
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I. INTRODUCTION 

A well-known result of BCS theory is the variational wavefunction, describing the ground state of a superconductor. 
In the limit Q — > this function can be easily extended to describe a superconductor [^j^ moving at a small and 
uniform velocity v = {2me)^^hQ 
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This function has the mathematical shape of a Bose condensate of pairs, where the wavefunction (j){f) describing the 
relative motion of electrons forming a pair is the Fourier transform of {[1 + (efe/Afc)^]^^/^ — (efc/Afe)}, and exp {iQ ■ R) 
is the macroscopic wavefunction describing the center of mass motion of each pair. The similarity to a Bose condensate 
wavefunction is somewhat misleading, as also the wavefunction of a gas of uncorrelated fermions can be written in 
this form, in which case (f){r) is a non-trivial function with an tail. In the limit of a weak effective interaction 
(j){r) has an algebraic tail just as for the free electron gas. If the interaction is strong, (j){r) can be interpreted as 
a wave-function describing the relative motion of two electrons forming a Bose-condensed pair. Q If the effective 
interaction is an on-site attraction, the electrons pair-up in a singlet-wave function with an enhanced probability to 
occupy the same site. Clearly if the electrons experience a strong on-site repulsion, the tendency towards pairing 
disappears. With a net attraction between electrons occupying neighbouring sites in the lattice, it is still possible to 
form a paired state, but (f)(r) has to be constructed such, that the particles avoid the same site. This condition is for 
example fulfilled when when (j){f) has a finite angular momentum. 

One may wonder whether the analogy to Bose-condensation can be drawn further, and consider the energy spectrum 
of pair-like excitations as a function of pair-momentum. This problem was first treated by Bogoliubov jj], and 
Anderson jsj. If the electrons experience an on-site repulsion, with a nearest-neighbour attraction, the collective 
mode spectrum becomes particularly rich. It turns out that several branches occur below the continuum of quasi- 
particle excitations, corresponding to order-parameter fluctuations of various symmetries of pairing The existance 
of low-lying collective modes may be important when attempting to identify a superconducting gap in the infrared, 
Raman, or inelastic neutron scattering spectra of these materials. 

Collective modes in superconductors have in the past attracted the attention for a variety of reasons: (1) Bogoliubov 
predicted the existance of a longitudinal collective mode with a sound- like dispersion B. Long range Coulomb 
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interactions make the spectrum identical to the plasmons of a normal Fermi gas, as was shown by Anderson 
(2) The collective mode spectrum naturally follows from a gauge invariant formulation of BCS theory 1^, and a 
consistent explanation of the Meissner effect requires that the whole interaction Hamiltonian (as opposed to the 
reduced BCS Hamilitonian) is taken into account [gj^]- (3) As collective modes mediate electron-electron interactions, 
plasmons ]^[ic|] and spin-fluctuations have been considered as possible candidates for a pairing-mechanism. 

(4) Certain modes, in particular condensate phase-fluctuations near or below the pair-breaking gap, are important for 
the thermal behaviour, notably Tc, of the superconductor [T^ ]. (5) An instability of the ground state and an incipient 
phase transition to a state with a lower energy follow from the softening of collective modes ||lq-fi8|. (6) Collective 
modes may show up in experimental spectra, such as in optical Jl9|-^ or Raman spectroscopy ]22| , ^ . (7) As there 
is no inter-plane hopping in a layered electron gas, the fc-dependent plasmon spectrum becomes gapless [ p4| , which 
may give rise to an interesting behaviour in the region for momentum and frequency values where the collective mode 
crosses 2 A p5| , p^ . (8) If there exists an electron-electron interaction in channels with a finite angular momentum L, 
excitons with the corresponding symmetries can exist P Jl^ , p^ . 

Usually modes of phase and density are treated separately from amplitude modes p^ , p3t , and spin-fluctuations 

|P^ , |ll|j26| . As we will see below, especially for a non- vanishing momentum a coupling exists between the four collective- 
mode channels of spin-density, charge, phase and amplitude of the order parameter. The aim of this study is, to derive 
general expressions for the collective modes in the superconducting state, using a unifled approach including effects 
of finite momentum pairing. In the last section examples are given for the collective modes and the generalized 
susceptibility in the superconducting state. It is shown that a d-wave superconductor may become unstable with 
respect to the formation of a spin-density wave, or possibly a mixed SDW-wave plus superconducting state, if an 
on-site repulsion is taken into account in addition to having an attractive interaction in the d-wave channel. More 
detailed calculations of various response functions and the comparison thereof to measurements on specific materials 
will follow in a future publication. 



II. THE MODEL HAMILTONIAN 



In the discussion of the collective modes we will make extensive use of two-particle creation operators. We will see 
below, that the channels with Sz = ~l, Sz = 1 and Sz = are decoupled. In the Sz = ±1 channels there are triplet 
pair-excitations, and spin-fluctuations. In the S'^ = channel there are spin-fluctuations, density fluctuations, and 
fluctuations of phase and amplitude of the order-parameter (singlet and triplet pair-excitations). The corresponding 
operators are in the same order 

(j)k{Q) = C-k-QiCkT - cl^Q^clf^^ 
ipkiQ) = c_fc_QxCfct +4+QTC-fei 

The remaining 6 combinations are c\_^Q^Ck(y^ (spin-fluctutuations with Sz — ±1), and C-k-QaCka with the corre- 
sponding hermitian conjugates (spin-triplet phase- and amplitude-fluctuations). When transformed to Euclidean-space 
representation these operators acquire a more transparant physical meaning. For example the spin density distribu- 
tion function n-^{r) — n^lr) has as its Fourier transform ak{Q)- Similar relations exist for the other operators, and 
the notation a{Q), p{Q), 4>{Q), and ip{Q), will be used to indicate the Fourier-transforms of the spin-density, charge 
density, phase and gap-amplitude distribution in Euclidean space. We will consider a system of interacting electrons 
which can be described with the following Hamiltonian: 

H = J2 Cfc/^fc(o) + E {\viQ)p(.Q)pi-Q) + luiQ) [p{Q)pi-Q) - ^(Q) ■ ^(-Q)]] (3) 

k Q ^ J 

where in V{Q) I lumped together the Coulomb interaction with all other spin-independent interactions, which could 
be due to the coupling of the electrons to the other degrees of freedom of the solid. In principle, and in particular if 
the interaction Kernel is derived from boson-exchange models such as electron-phonon interaction, there can also be a 
separate dependence on the momentum of the interacting particles. For compactness of notation I will not explicitly 
include such a k and q dependence in the Hamiltonian. 

With the spin-dependent interaction assumed here, the total spin of the system is still a good quantum number. Such 
terms can appear if the model Hamiltonian is derived from a more fundamental one by projecting out part of the 
Hilbert space. A well-known example is the occurance of the Kondo exchange-interaction in a magnetic impurity 
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system after carrying out the SchriefFer- Wolff transformation. Other examples where such terms occur are the RKKY 
interaction in magnetic alloys, and the superexchange in rare earth-doped semiconductors. Also the on-site Hubbard 
U term is usually written in this form, although in this case the Pauli-principle already automatically excludes occu- 
pation of the same site with parallel spins. 

As we will discuss the equations of motion of the collective modes for a general form of the effective electron-electron 
interaction, it is worthwhile to summarize the expressions for the gap equation and the free energy. The thermody- 
namic potential at T = of a BCS superconductor is the expectation value of the grand canonical Hamiltonian, and 
is easily obtained by taking the expectation value of Eq. ^ using the variational wave function of Eq.|l| 

V, Vk,,...,Vk^) = 2J2H\^ (Cfe -^^)+Yl {ukVkXkqU*v* + \uk\^Vik - q)\vg\^) (4) 

k kq 

where Xkq is the pairing potential. For the type of interaction introduced above one obtains Xkq = V(k — q) + ^U{k — 
q) + ^U(k + q) where the last term is the spin- flip scattering contribution contained in (Jk{Q) ■ <^q{^Q)- The last term 
in Eq. corresponds to the exchange energy. From f2 one obtains the gap-equation by calculating the minimum as a 
function of the set of variational parameters Vk^, ...,Vkpi- The number of particles in the ground state is obtained by 
taking the first derivative with respect to fi. The resulting set of equations is 



l^a 'J-q'Jq^kq — |„.|2_|^ 



Apart from a shift in chemical potential the effect of the exchange energy term on the thermodynamic potential is to 
renormalize the single particle dispersion Sk ~ ^k ~ M; which now has to be replaced with ik = ^k — l^qP^C^ ~ 
After the ground state has been found from minimalization of the free energy, the quasi-particle spectrum is obtained, 

with the usual BCS-type energy dispersion Ek = [e^ + A|] ^^'^ , and with Ak/Ek defined as 2ukVk- In the following 
sections I will also use the (standard) notations bk = UkVk and Zk = — \vk\'^)/2. 

If Xkq has a non-trivial fc-dependence we can make a partial-wave decomposition 

^kq = ^llJa{k)KllJa{q) 

a 

where {ipaik)} is a complete set of orthogonal functions, chosen such as to diagonalize the pairing potential. We can 
make a similar expansion of the order parameter Aj, — A^^paik) with the help of which one obtains the coupled 
gap equations 
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We notice, that for A ^ a decoupling of pairing-channels occurs, depending on the presence of off-diagonal elements 
in the decomposition of l/E'fe ^ l/|efc|. As A^^ is real, the set {-00; (A;)} can be chosen as real numbers. As a result 
also Aa is real. Solutions like "s-l-id" |^ become possible if there is a degeneracy between solutions with a different 
symmetry. 



III. EQUATIONS OF MOTION 

The equations of motion are of the form [H, O] = vO, where O is a linear combination of pair-operators representing 
an excitation of the system with energy v. Although these equations have been treated extensively before, in the 
previous papers the coupling to the collective spin oscillation channel has not been considered. In particular a spin- 
dependent term was not included in earlier publications. As one of the aims of this paper is to discuss collective modes 
of spin-density in the superconducting state, I re-derive the equations of motion with this extended hamiltonian. 
In the superconducting state the equations of motion of spin density((jfe(Q)), charge density (pkiQ)), order parameter 
phase {(pkiQ)), and order parameter amplitude (ipkiQ)) are coupled in a non-trivial way. The commutator of each 
of these two-particle operators with the interaction part of the Hamiltonian generates products of four single-particle 
operators, which are approximated by taking the expectation value of all combinations of two of the operators 
appearing in this product. The resulting terms fall in two categories: those which have the same fc-value, and those 
which are a weighted summation over fc-space. The latter give rise to the collective modes. In the first category one 
obtains (1) self energy terms which can be absorbed in a shift of the chemical potential, (2) exchange self energy 
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terms, due to which Ck is renormahzed to ik = J2q I'^gP^C^ ~ ^'^^ (3) cross-terms proportional to Afe, hnking 
CTfe to ipk, and pk to ^^-operators. 

Finally the category of weighted averages of two-particle operators over fc-space involves both direct and exchange 
terms, and is given by the expressions 



where I introduced 



Sfe(Q) = EgHUk^q.QKiQ) 

Rfe(Q) = EgHj,{k,q,Q)pgiQ) 

AkiQ) = EgHlik,q,Q)cj,,iQ) 

Bfc(Q) = EgH'i,{k,q,Q)MQ) 



Hl{k,q,Q) = -\U{Q)-V{k-q)-\U{k-q) 

Hl{k,q,Q) = 2V{Q) + h^U{Q)-V{k-q) + ^U{k-q) 

Hlik,q,Q) = V{k~q) + \U{k-q) + \U{k + q) 

^{k,q,Q) = V{k-q) + \U{k^q) + \U{k + q) 



(6) 



(7) 



With these definitions, and using the random phase approximation described above, the commutators of the pair 
operators can now be derived. The actual calculation is a straightforward, though rather laborious, exercise in 
commutator algebra. A detailed description of the various terms has been given by Anderson, and later discussed 
more extensively by Bardasis and Schricffcr, who retained a number of vertices in their final analysis which were 
neglected by Anderson. In the present paper all vertices discussed by Bardasis and Schricffcr are taken into account. 
The expressions are however modified due to the spin-dependent interaction term in Eq. ^. The set of commutators, 
including the exchange interactions, is 

[H,ak{Q)] = -e^QPkiQ) - ^IqMQ) +z^Qnk[Q)-blQBk{Q) 

[H,pkm = IkQ^kiQ) - KqMQ) +z^QSk{Q)~b+QAk{Q) 

[H,MQ)] = -~4qMQ) - ^tqPkiQ) -btQnk{Q)^ztQ^k{Q) ' 

[H,MQ)] = ~~4qMQ) ~ ^kQ<^k{Q) -b-^SkiQ) ~ z+QAk{Q) 

^k,£k, bk and Zk were already defined in the previous section. For the sake of compactness of notation I introduced 
^kQ = bk+Q ± bk, = Zk+Q ± Zk-, A^g = Afe+Q ± Afc and e^g = ik+Q ± ffc- 

The first two terms of all four commutators correspond to (1) the kinetic energy with exchange self-energy corrections 
(Fig. |l|a), and (2) Boguliobov-Valatin particle- hole mixing (Fig. |l|a'). The remaining two terms in each of these 
expressions can be better described with reference to the definition of the collective coordinates in Eqs. ^ and 0. 
Let us first consider Ra;(Q) and Sk{Q)- The V{Q), U{Q) and U{k — g)-terms correspond to the polarization vertex 
in the commutators of ak and pk (Fig. |^). In the commutators of 4>k and ipk the V{Q), U(Q) and U{k — q)-term is 
a polarization vertex combined with a particle- hole transformation on one of the legs (Fig. l|b'). The V{k — q)-teTnia 
correspond to the exchange scattering vertex without (commutators of Cfe and pk,Fig. He) and with particle- hole 
transformation (commutators of (pk and V'fc, Fig. ^j'). 

Finally Ak{Q) and 'BkiQ) correspond to the direct particle-particle scattering vertex without (commutators of (pkiQ) 
and ipk{Q), Fig. |l|d) and with particle-hole conversion (commutators of <Tk{Q) and pk{Q), Fig. 0d'). 
If we apply the equations of motion to a general operator of the form 

= Y1 Kk{Q)'^k{Q) + v2M{Q)Pk{Q) + v3.^k{Q)MQ) + Vi^k{Q)MQ)] 

k 

we find that they can be written in matrix form as 

H\k,Q)vk{k,Q) + Y,H\k,q,Q)T{q,Q)vg{Q) ^ vvkiQ) (9) 

9 

The interaction Hamiltonian iJ' contains the matrix elements of Eq. ^ on the diagonal, and is zero elsewhere. We 
furthermore use the zero'th order Hamiltonian describing non-interacting quasi-particles 
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(10) 



4 



and the dimensionless matrix containing coherence factors 



T{k,Q) = 



( ' 


^kQ 










-Kq 


V ~Kq 






-Kq \ 
-Kq 
-Kq 



(11) 
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The collective modes can be found by looking for poles in the correlation functions, in particular the density-density 
and the spin-spin correlation functions << T p(f, t) p(r' , 0) >>,y and << Ta{r,T)a{r' ,0) >>u, where p{r,T) etc. are 
the Heisenberg representation of the operators exp iQ ■ rp{r) which were defined in the previous section. In the 

superconducting state also << T(j){r,T)(l){r',0) >>i, and << Ttp{f,T)4>{r' ,0) »i, become relevant. Together with 
the six off-diagonal correlation functions, the four diagonal functions define a 4 x 4 two-particle Green's function 
matrix. The matrix Kl^{Q^v) = {v — H^{k,Q) — iQ^)^^5k,q corresponds to the Lehmann representation of this 
Green's function in the absence of residual interactions (i.e. with W = 0). As this describes the response of a gas 
of non-interacting quasi-particles there are no poles corresponding to collective modes. The generalized susceptibility 
X°(Q,i') = J2k q^i^^Q)^k qiQ^'^)- 1^ ^'-'^ include H^, we can calculate the Green's functions in the RPA by 
applying the Dyson equation 

KkqiQ, v) = Kl^{Q, v)5k,q + Klk{Q, v) J2 H\k, k', Q)T{k', Q)Kk'q{Q, v) (12) 

We can use the same partial wave decomposition as introduced in the previous paragraph where we discussed the gap 
equation. It is straightforward to show, that the above Dyson equation has the solution 

Xc.AQ^ ^)^Y. ^) (1 - hW^I (13) 
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where I used the partial wave decomposition 

Xa,/3(Q, i^) = X! i^a{k)T{k, Q)Kkq{Q, J^)'0/3(g) 
kq 

with similar expressions for x^, and H^. The collective modes correspond to the zero's of the determinant of 

So.jS,,j (14) 

which can be determined numerically, and in some limiting cases also analytically. The expression of the response 
function Eq. |l^ corresponds to calculating the series of diagrams depicted in Fig. |l|. It is possible to improve further 
by taking into account the screening of the vertex in all of these diagrams, except in the polarization vertices of 
Figs. |l]b and |l|b', as this would lead to double counting of the vertex corrections. (N.B.: Although in this paper the 
pairing interaction is introduced as an independent model parameter, one should keep in mind, that for an electronic 
mechanism of superconductivity such as a spin-fluctuation or plasmon-intermediated interaction, the pairing arises 
precisely from these vertex corrections.) This procedure was proposed by Anderson, Rickayzen and also by Bardasis 
and Schrieffer. Moreover, in the next section we will see, that in the normal state the a and p channels are completely 
decoupled for all values of Q. This implies, that the sum over diagrams for the charge fluctuations does not contain 
any vertex correction due to the spin fluctuations and vice-versa. Hence, it is necessary in this case to screen all 
vertices in the charge-fluctuation channel with the spin- fluctuations, and vice versa. As has been shown by Rickayzen, 
in the superconducting state the screening properties are basically the same as in the normal state 0|. 
One has to be cautious with this procedure of screening the vertices, as, by making the RPA before calculating the 
sum over diagrams, certain classes of vertex corrections are omitted. As a result inconsistencies may arise, as can 
be seen from the following example: If we consider the Hubbard U model, the on-site interaction can be introduced 
either using an on-site spin-independent (V) or a singlet-only (U) term as defined in Eq. ^ The expressions for the 
equation of motion should be independent of this choice, as the Pauli-exclusion principle automatically projects out 
the double occupancy of the same site with equal spins. Indeed, we can check from Eq. ^ that this requirement is 
satisfied as long as we do not introduce screening. If we follow the recipe, that in the first two lines of Eq. |^ the 
polarization diagrams U{Q), V{Q), and U{k — q), but not the exchange diagram V{k — q), should be replaced with 
the bare interaction, we arrive at a different result depending on whether we introduce the on-site interaction through 
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a singlct-only or a spin-independent interaction. 

This inconsistency is removed, if we replace the direct and exchange terms in with the charge screened value. In 
the same way screening with spin fluctuations should be introduced 'by hand' in the direct and exchange terms in 
Hp. Finally all three terms in and H^^ should be replaced with the charge and spin fluctuation- screened vertices. 
Let us now calculate K'^ by inverting [i/ — H'^p] . The determinant is 



1^ _ = ^4 _ 2^2(^2^^ + ^2) + (^2^^ _ ^2^2 



(15) 



The zero'th order two-particle Green's function is then 



-f+A + A- 
-i/e-A+ 
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+ (?+?- + A+A-)A- 

i/e+A+ 
-i/e~ A~ 



-vt-A+ 
-i/=A+ 

+ (e+e- + A+A-)A- 
+e-A+A- 



e+e-A+ 
+(A+^ - i.2)A- 



i/e+A+ 
-i/e~ A" 



+e-A+A- 



(16) 



The 4x4 matrix K^T becomes 

K°T = \u - H""]-^ 
I -v^{z-ir + b^A^) 
+ (&-A+ + 0-e+) 
x(e+e- + A+A") 



+!/«-(?+ +A- ) 



z V 

+vz~{e+^ + A+^) 
+vb+(e+A- - e-A+) 

-iy2(z-£- +6+A+) 
+ (b+A- + 



iy2(b+f- - Z + A-) 
+ (z+A+ - b+e+) 
x(e+e- + A+A") 

b+i/3 



b-!.3 

+i/0+(e+A- - e-A+) 
-ub-{e+^ + A+^) 

;/2(b-J- - z+A^ 



+vb-{€+A+ - e-A-) x(e+e-+A+A-) 



+i/«+(e+A+ - e-A") +(z+A- - b-?+) 



iy2(-b"f+ + Z-A+) 
+ (2-A- +b-e-) 
x(e+e- + A+A") 

-i/2-(e+A+ - e-A") 
\ -1.6- (e-^ + A+^) 



-1.6+ (e+ +A- ) 

-(y2(2 + £+ + b+A+) 
-i.2-(e+A- - e-A+) +(6+A-+z+e-) 
_i.6+(e-% A-^) x(e+e-+A+A-) 



1.^6+ 



-i.2(6+e+ - Z-A-) 
-(«-A+ - b+e-) 
x(e+e- + A+A-) 



-1.2+ (e- +A+ ) 
+i.6+(e+A+ - e-A-) 



<(?+?- + A+A-) 

1.32 + 

-i.2+(e-^ + A-^) 
-i.6-(e-A+ - e+A-) 

-u^{z+e+ + b'A-) 

+ (2+e-+6-A+) , 

x(e+e-+A+A-) / 



(17) 



From inspection of the matrix elements it turns out, that they all contain the factor (z/^ — {Ek+Q — Ekf') in the 
numerator. As the same term appears in the denomenator, these factors cancel. K^T turns out to be symmetric, and 

the exact result is 



TK" = K"T = 



I —Ek+qEk + A^Ak+Q 



^, ^k+Q'^k-^k^k+Q 
'^k + Q+'^k 

^fc+gAfc — EkAk+Q 



"-^k + Q^k-^k^k + Q 
'^k + Q+'^k 



'^k + Q^k-'^k^k + Q 
Bk + Q+l^k 

-^k+Q^^k — ^k^k + Q 

+^kf-k + Q 

Bk + Q+'l'k 



-Sk+Q^k — ekAk+Q 



Sk + Q^k — Sk^k + Q 



. ^k+Qf^k + ^k^k + Q 



-Ek + QEk - Afc + gAfe 
— tk + Qf-k 

^k + Q^k + ^k^k+Q 
^k + Q + l^k 



^k + Q^k-^k^k + Q \ 
^k+Q + l^k » 



-efe+gAfc — efeAji^Q 

^. ^k + Q^k+^k^k + Q 
^k + Q+^k 



-Ek+QEk + AfcAfc+Q 



(18) 



IV. EXAMPLES 



In this section I will apply the formalism outlined above to a number of examples with an increasing degree of 
complexity vis a vis the type of electron-electron interaction that is assumed. The energy dispersion is assumed to be 
of the form 



6 



s: 


Mk) 


= 1 


s*: 


Mk) 


= cos(fcj:a) + cos{kyb) 




Mk) 


= cos{kxa) — cos{kyb) 


Px- 


Mk) 


= ^/2sin{kxa) 


Pv- 


Mk) 


= \/2sin{kyb) 


dxy- 


Mk) 


= 2sin(kxa)sm{kyb) 


etc. 







efe = — 2t {cos{kxa) + cos(fcj,a)) — 2t' cos{kxa) ■ cos(fcj,a) — fi (19) 

where a b and c are the lattice parameters. The t and <'-tcrms are due to nearest-neighbour and next-nearest neighbour 
hopping in a square lattice. If t' = at half filling of the band, such a dispersion relation has the remarkable property 
that the Fermi surface forms a perfect square, with a diverging effective mass over the entire Fermi surface. In practice 
this situation will never occur, as there will always be some finite coupling between next nearest neighbours. This 
causes a bulging of the Fermi surface, which eventually transforms into a rotated Fermi surface if \t'\ ^ \t\. 
In all examples I will restrict the discussion to systems where electrons have an on-site attraction or repulsion, a 
nearest-neighbour interaction, or both, as well as the long-range e^/r repulsive interaction. Moreover the discussion 
is limited to the situation where a single band crosses the Fermi surface, and tight-banding language will be used for 
the description of this band. In particular I will consider a tightbinding band on a three-dimensional square lattice, 
with a strong anisotropy leading to quasi two-dimensional behaviour. A convenient set of functions to be used for the 
partial- wave decomposition of H'' is then the set of harmonic functions: 



(20) 



The /c-space representation of the on-site Hubbard U interaction UoUi^riii is the fc-independent function U{k — q) = 
Uoipo{k)'ipo{q) ■ If we consider the nearest-neigbour interaction ^Vi X]<i j> "-i^j we find, by means of a direct Fourier 

transformation of the operators rii = cl^Cii + cl^Cii, that this can be cast in the form ^J2Q^iQ)piQ)pi~Q) with 
V{Q) = 2Vitpi{Q), so that we obtain the partial-wave decomposition 

Vik -q)^V, iMk)MQ) + Mk)Mq) - Mk)Mq) - Mk)Mq)) (21) 

A singlet-only nearest-neigbour interaction jUiJ2<ci j>i'4^'^'ji ^ ^li'^l^)i'^jl^^i1 ^ •^jT'-u) be cast in the form 
I J2q U{Q) [p{Q)p{—Q) — <y{Q) ■ (?(— Q)] with U{Q) — 2Ui'4'i{Q), hence it has the same partial wave expansion as 
Wiipi^Q). However, from Eq. |^ we see, that the singlet-only interaction has other prefactors, and is summed over 
U{k — q) and U{k + q) in the pairing channel. 

Finally we have to take into account the long range Coulomb interaction. Here we will use the lattice Fourier trans- 
form of e^/r. The screening of the Coulomb interaction for part of the vertices has been discussed above, and is 
essential, as a bare interaction is known to create a singularity at the Fermi level within the random phase 
approximation. I will use the convention in the remainder of this paper, that V{Q) is the bare Coulomb repulsion at 
large distances, whereas for shorter distances Uq and Vi are the projections of V{Q) on the on-site interaction and the 
spin-independent nearest neighbour interaction respectively. Taking all these terms together we obtain for a model 
with a 'singlet-only' nearest neighbour interaction 

ij;(o,o) = -u^-u^MQ) 

Hl{a,a) = -C/f/2 (« = 1,2) 

i7;(o,o) = 2v-{Q)~us + u^MQ) 

Hl{a,a) = C/f/2 (« = 1,2) ^^^> 

ij;(0,0) = U^^ 

ij;(a,a) = (a = 1,2) 

For all symmetries we have H^{a, (3) = H^^(a, [3). The upper-indices p and a indicate whether screening with charge- 
or spin-fluctuations is implied. The minus-sign in front of the Uq term in i/*(0, 0) is not a misprint. As 2V{Q)^ 
'contains' the on-site Hubbard term, the sum of these two contributions is +?7o. In principle one should also include 
higher harmonics, as the expansion of V{Q) does not end at V'4- However, as the expansion only appears as a screened 
interaction in expression it is reasonable to work with a a model where such interaction-terms are neglected. If 
the nearest-neigbour interaction is spin-independent we must also include px and py symmetries of pairing, and we 
obtain 
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hUo,o) 


= -u^ 

u 








= -vr 


{a 


= 1..4) 


ij;(o,o) 




-US 






= -vr 


{a 


= 1..4) 


^^;(o,o) 


— Uq 








= 


{a 


= 1..4) 



(23) 



In the previous section we have seen, that in addition to the partial- wave expansion of i/', we also have to make a 
similar expansion of TK^. The expression for this product is given in Eq. |l8|. The partial- wave expansion of this 
expression is in general complicated, and has to be done with the help of a computer. Some limiting cases exist 
however where the integrals can be solved, especially when an expansion for small Q can be made. Some of these 
limiting cases will be treated in the subsequent sections. In addition numerical calculations will be given at general 
values of the collective mode momentum Q. 



A. Normal state limit 

In the non-superconducting limit Eq. ^ has only non- vanishing matrix elements on the diagonal. Furthermore only 
the charge and spin-channels are relevant in the absence of off-diagonal order. Let us make the further assumption 
that the electrons interact with each other via an on-site Hubbard U repulsion, which is therefore independent of k. 
After the summation over k we obtain for the top-left corner of Eqs. ^ 

where the fk are Fermi occupation factors. Let us assume for this part of the discussion, that the mo- 
mentum of the electrons in the plane is unbounded. In that case the Fourier transform of e^/r is dis- 
crete in the direction perpendicular to the planes, and continuous along the planes, so that |Q V{Q) — 

27re^(i^ sinh (Qyii) ^ [cosh ((5||d) — cos (Q^d)] ^ where d is the interlayer distance. Let us define Q = 

Aire^ /V{Q), which has the property limg^o Q — Q- If we assume that we have a cylindrical Fermi surface, with an 
isotropic Fermi velocity vp, and a Fermi wavevector kp, we obtain with 1 — e = 2F(Q)x"(2, 2), and I'p = 2e^d~^hkFVF 

-1/2 ^ (25) 

The plasma dispersion relation becomes 



I /-, H'QivF 



The spin-susceptibility per unit cell (f2„ is the area in the 2D plane) is 



(1 - n2Q|z;|/i.2)-l/2 _ 1 
W + Ua[{l^ri^Qlvl/v^)-y^ - 1] 



where W = is the effective bandwidth. We see, that in the high frequency limit [v 3> Tt-QwVf) X^^(1,1) = 

i?FQ|f^M[i^^ + UoEpQ'^iyiu]^^ ^ and in the low frequency hmit x^-'^(l, 1) = {Uq — W)^^, Hence the AC susceptibility 
is suppressed, whereas the static susceptibility is enhanced. A magnetic instability occurs for Uq ~ W . The above 



expressions are derived assuming a free electron dispersion. If the Fermi surface has nesting vectors |12|, instabilities 
for specific values of Q are often found. 
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B. s-wave superconductivity 



For s-symmetry, and neglecting the radial fc-dependence of the pairing potential, the partial wave decomposition of 

x'^ and is trivially achieved by summing over all k. As a model for the pairing-interaction wc adopt U{k — q) = —g. 
As is ususally done in the gap-equation, one can limit the energy-range of the interactions in these expressions by 
putting A = for energies larger then a scaling value (the Dcbije frequency for phonon-mediated pairing). For the 
long range Coulomb interaction we take again V{Q) = Ane^Q^^. Due to the fact, that e/c+Q = £-k-Q we obtain after 
summation that x°(l)2) = x°(l,3) = x°(l)4) = 0. Hence the spin-fluctuations are fully decoupled from the other 
three and can be considered separately. The remaining diagonal and off-diagonal susceptibilitues are finite, and the 
following expressions are required 

■Y^f-l 1 "1 = (Ek+Ek+Q){EkEk+Q—AkAk^Q—ekek+Q) 
2EkEk+Q{i^^-[Ek+Ek+Qf) 
O/o 9") = {Ek+Ek+Q)(EkEk + Q+AkAk + Q — tk€h+Q) 

^ ^ ' > 2EkEk+Q{u^-[Ek + Ek+Q?) 

O/o Q"l = {Ek+Ek+Q)(EkEk+Q+A.kAk+Q+ek^k+Q) 
^ ' 2EkEk+Q{u^-[Ek+Ek+Qf) 
-yO(A 4") = {Ek+Ek-^Q)(EkEk + Q-AkAk+Q+tkik-\-Q) 

^ ^ ' > 2EkEk+Q{y^-[Ek + Ek+Q]'^) 

a — (Ek+Ek+q) C28) 

2EkEk + Q{u''-[Ek+Ek + Q?) 
rp y-v {Ek + Ek + Q){ck-\-f..h^Q) 

~ ^k 2EkEk+Q{y^-[Ek+Ek+Q]^) 
_ (-Efc-l--Efc+o)(efc-efc+Q)^ 
^k 2EkEk+Q{u^-[Ek+Ek+Qf) 

^k 2EkEk+Q{i'^-[Ek+Ek+Qf) 



Let us first consider the spin susceptibility. As now x(l) 1) = i-gj^'(ii) ^it^ 5 > for a BCS interaction, the 
spin-fluctuations are pushed to a slightly higher energy. 

The generalized susceptibility for density, phase and amplitude can be expressed using the above definitions as 

/X°(1,1) + 2A25 -lyAS AT \ 

-uAS x°(3,3) -i^{T + N)/2 

\ AT -y{T + N)/2 x°(4,4) / 

Using the gap equation, (1 = Y.k ife) "^'^^^ *° prove, that x"(3,3) = -l/g + v'^S/2 - M/2 and X°(4,4) = 
X°(3,3) — 2A^S. Using these properties, we see that the matrix 1 — -ff*x° becomes 

/ 1-2Vqx°(1,1)-4FA25 2VqvAS -2VqAT \ 

-gvAS g{v'^S-M)/2 -gu{T + N)/2 

V gAT -gv{T + N)/2 g{y-AA'']S-M)/2) 

and the collective modes can be calculated from the determinant 

= (1 - 2t/QX°(l, 1)) {{v^S - M) _ 4A2]5 - M) - i,\T - Nf) 

+4yA2 {{N^ + 2TN)Sv'' + M{S''y - 4A2] - T^) ^^^> 

To further analyze this expression we need to make a series expansion for small Q. In what follows we will neglect 
A'' (oc Q^, but with a vanishing prefactor if the gap has electron hole symmetry). Furthermore we notice, that we 
can write 2T/qx°(1, 1) « VpiQfiv'^ -W^ - 4A2)-i, Af « 5 < (uf • Q? >, VqM w (1 + ly V4A2)Vqx°(1, 1), and 
T « —2\xfS, where [If = Ep — W/2, and W is the effective bandwidth. Retaining only leading orders in < {vp ■ Q)"^ > 
we see, that the collective modes can now be solved from 

2(1.2 _^r2_ 4^2 _ j,^(g)2)^j,2 ^ 2A2) < {vF ' Qf >= iv^ " 4A2 - 4/x|)(i/' -W'' - 4A2)(i/2 - v^{Qf) (30) 
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Hence we see, that the negative U Hubbard model permits in principle four collective modes: a spin-density oscillation 
discussed above, a plasma-mode, and two additional modes, which are however situated in the two-quasiparticle 
continuum, and therefore are strongly damped. Interestingly the plasma-like mode can exist at frequencies above 
and below the gap, depending on the initial value of Vp{Q) in the normal state. As has been discussed by Fertig and 
das Sarma pst a layer dispersion relation as discussed above, permits the existance of low lying plasmons below the 
gap. Another mechanism for reducing the plasma frequency in the superconducting state is strong damping of the 
motion perpendicular to the planes, as we recently discussed E^. 



C. s-wave superconductivity in a layered electron gas. 

If the superconductor is strongly anisotropic, the plasma energy for Q ^ depends on the direction of propagation. 
An extreme example of this arises, when the mass in one of the three directions is infinite, resulting in a system which 
behaves two-dimensional from the point of view of the signle particle band-structure, whereas the Coulomb forces are 
three dimensional. A simple model exhibiting such behaviour is an infinite stack of a two-dimensional layers. The 

using hydrodynamic calculations, as 
well as with the random phase approximation. The resulting plasmon spectrum of such a metal is, in the limit of a 
large wavelength, v{Q) — VpQ\\\Q\~'^ , which for Q±^ = saturates at the value Vp, while for finite values of Q±^ it has 
an accoustic-like dependence on Q\\. 

This implies that here we have a system which on the one hand has a density of charge carriers characteristic of a metal, 
and, provided that there is a pairing mechanism, therefore has the potential of becoming a BCS-like superconductor. 
On the other hand the dynamical response of the electrons in one of the directions is more characteristic of a 
semiconductor or an insulator. This combination provides us with an example where the Anderson-Higgs mechanism 
does not shift the Goldstone mode to a high energy, in spite of the fact that the particles interact through a long 
range Coulomb force. Here we will use the dispersion introduced in Eq. |l^ with t' — 0. In this example W — At is the 
bandwidth. For the long range Coulomb forces we take the lattice Fourier transform of e^/r, which has ^~^A'neQ~^ 
as its long wavelength limiting behaviour, and posesses the same periodicity in fc-space as the tight-binding band. 
Ry* = e^d~^ is the effective Rydberg which, together with the Fermi energy, sets the scale of the plasma frequency 
in the planar direction [vp = 2a/ EpYly* for Ep « W). For the pairing interaction we adopt the same model as in 
the previous section, with the following set of parameters: Ep/W = 0.35, g/W = 0.6729 (resulting in A/W = 0.25), 
and Ky* /W = 4.0. In Fig. ^ the result is displayed with Qj^ as a parameter in the range from to d^^n, and 
agrees well with the calculations of Fertig and das Sarma, and those by Cote and Griffin ||2^,|l^ . Due to the model 
assumption of an energy independent attraction, the electron-hole continuum becomes a broad band already for zero 
momentum in the two-particle channel. In Fig. ||b the same set of calculations is displayed for the normal state. 
The plasma frequencies become somewhat smaller in the superconducting state, which is due to the fact that the 
gap in this example is relatively large. It reflects a well-known property of the negative U Hubbard model, that 
the mass of a pair is enhanced due to the fact that two particles have to hop simultaneously Also a strong 

qualitative difference arises, which is not directly evident from these curves. This is the change in character of the 
modes. In Fig. ^ the distribution of weight of the mode over the density, gap-phase, and gap-amplitude branches 
is displayed. First of all we notice that the contribution of the latter is neglegible. The second interesting feature 
is, that the nature of the collective mode changes gradually from a pure phase-fluctuation at low energy to a 50/50 
percent phase-density mixture at the edge of the particle-hole continuum. Inside the electron-hole continuum the 
collective modes are damped. (Although they may still persist as a resonance, they can not be indentified from the 
zero's of a determinant) However, for energies larger than the particle hole continuum of our band we see that the 
density-fluctuation character dominates. 

The plasmon-dispersion, which is coupled to the s-phase fluctuating channnels for small q, is accoustic for finite 
Q±/Q\\ so the Landau criterion [ p8| remains satisfied in spite of having a gapless plasmon spectrum. 



electrodynamics of this system was already discussed by several authors 1 24 



D. From s- to d-wave superconductivity in a layered electron gas. 

Let us consider a spin-dependent nearest-neighbour pairing interaction (Ui). As a result the only channels open 
for pairing are s and d-wave, and these are also the only channels in which a bound state can be pulled below the 
electron-hole continuum. Let us furthermore assume the same bandstructure as in the previous section. The phase- 
diagram, displayed in Fig. ^ was calculated by searching numerically for the minumum of the Helmholz free energy 
(Eq.^ -|- fiNe, using /i as a Lagrange parameter to keep the number of electrons fixed) as a function of A^, A^* and 
Ad(2;2_j,2'). The boundaries, which are indicated in this diagram are calculated for T — O.OIW, where W — At is the 
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bandwidth. For T = no sharp phase boundaries exist. Somewhat surprisingly, for \Ui\ larger than a critical value 
(which depends on rie), the ground state is of mixed as + j3d symmetry, which is automatically a state of broken 
spatial symmetry. However the strong interaction which is required probably does not exist in any realistic model of 
superconductivity. It is worthwhile to mention in this context, that the region of sd-mixing almost coincides with the 
region of p- wave symmetry, if we use a spin- independent interaction [Vi) instead. 

Let us now consider the collective modes for the examples along the line AB indicated in the phase diagram. The 
result is displayed in Fig. ^. We see, that a soft-mode developes if we approach the phase boundary between s 
and d-wave superconductivity. The transition takes place exactly when the mode has developed into a Bogoliubov 
sound-mode. If we keep imposing the s-symmetry for the ground-state, while actually being in the d-wave part of 
the phase-diagram, we always find a soft-mode of phase-fluctuating character, indicating that the solution is instable. 
If we allow the groundstate wave-function to become d-wave paired, the gap disappears, and a sound-wave phase- 
fluctuation mode occurs directly below the particle-hole continuum. Both the Bogoliubov mode and the lower bound 
of the particle-hole continuum are sound-like, so that according to Landau's argument [ p8| a supercurrcnt-flow is still 
possible in spite of the fact that there is no gap. 



E. Phase versus spin-fluctuating modes in a layered electron gas. 



Let us now consider the singlet-only nearest- neighbour pairing interaction Ui. In the discussion of the resonating 
valence bond state p9|-|3l[| the t-J model has been used, where J = —Ui, and a reduction of the double occupancy 
of the same site is included, either by replacing the bare hopping parameter t with an effective one, or by using 
more elaborate schemes. It is not the aim of the present discussion to address the tJ model. Instead we consider a 
Fermi-liquid, with an on-site repulsion {Uq) which is not too strong, and an attractive interaction between electrons 
on a neighbouring site {Ui). As the actual bandstructure in these systems is experimentally known to be better 
described by the three band model of Zaanen, Sawatzky and Allen (which is again a simplified version of the 
real valence band structure involving 6 oxygen 2p bands and 5 cupper 3d-bands for the occupied states, as well 
as unoccupied 3s and 3p states) a transformation to a single band hamiltonian will in principle generate both an 
effective Hubbard C/q and an intersite Ui |^ ^. Examples of such transformations can be found in the work by 
Emery [ p6| , and by Jansen p7[ |. However, also other, more complicated types of interactions are generated when 
making such transformations, notably the correlated hopping term (with six operators) which, as has been shown 
by Hirsch, promotes superconductivity of hole-carriers p^ . The interaction considered by Jansen as well as the 
correlated hopping term treated by Hirsch, effectively provide an on-site attraction, which, when considered on its 
own, promotes pairing in the (non-extended) s-wave channel. Also Ui term contains contributions from the virtual 
exchange of spin- fluctuations ||3|,|l^ . As has been discussed by Scalapino ||l2| , such processes give rise to an attraction 
on nearest neighbour sites, and increase the on-site repulsion between electrons. As the exchange spin-fluctuations are 
really vertex corrections due to the iJ^-channel, one could schematically regard t/i in Eq. ^ as the vertex correction 
of C/q. As such corrections are necessarily retarded, and therefore rather ill-represented by the non-retared interaction 
assumed here, the present analysis can at best provide a qualitative picture. 

BZA 30 1 considered pairing of the s*-type near half filling, Emery considered (ia;2_j^2 -pairing, and Kotliar studied 
both s* and d-type pairing. As we will see, the s*-type pairing is not a stable solution near half filling, and is 
dominated by pairing of the d-type. As the latter again tends to be unstable with respect to the anti-ferromagnetic 
Mott-Hubbard insulating state at half filling, superconductivity can only exist sufficiently far away from this region. 
As the optimal Tc would have been reached at half filling for a symmetrical band, this would lead to the conclusion 
that superconductivity is only a marginal effect in such a system. However, the high Tc cuprates do not have an 
electron-hole symmetrical band, and the Fermi surface is known to be strongly distorted from the perfect square that 
arises from considering only nearest neighbour hopping. This actually comes to rescue: As a function of band-filling 
it pulls apart the regions, where anti-ferromagnctism and high Tc have their highest stablility. 
The three coupled gap equations are (with x = k^a and y = kyo) 



E 



ta.TL\\Ek/{2kBT) 
2Ek 



Uq Uq [cos X + cos y] Uq [cos x — cos y] 

U 1 [cos X + cos y] Ui [cos x + cos y] ^ Ui [cos^ x — cos^ y] 
C/i [cos a; — cos y] C/i [cos^ a; — cos^ y] C/i [cos x — cos ^ 




(31) 



together with a fourth expression, which determines the chemical potential by constraining the electron-occupation 
number J2ki^ ~ ^k/Ek) = Ng. For a sufficiently small value of Ui or for T near Tc, where A becomes small, the 
denomenator has the four-fold symmetry of the crystal, and the cross-terms linking s* to d are zero for symmetry 
reasons. Hence only s and s* are coupled provided that Uq ^0. If A becomes large compared to the bandwidth, a 
priori there is no reason why mixing between s and d is forbidden, and indeed we will see, that such a mixing takes 
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place for a large value of Ui. 

I still need to specify the electron dispersion relation before we can solve the gap equations. For the dispersion relation 
we now use Eq. |l^ with t' — —0.7t. The shape of the Fermi surface obtained with this choice of parameters is very 
close to what has been calculated with the local density approximation for e.g. La2Cu04 and YBa2Cu307 po| , ^l| . 
Due to the finite value of t' a significant change occurs in the density of states (DOS) at the Fermi energy as a function 
of the number of electrons per unit cell. The DOS is now a-symmetric, and the maximum is shifted to the 'hole-doped' 
side of the point where the band is half filled. Of course the direction in which this occurs is dictated by the sign of 
t'. With t' < we mimic the situation encountered in the Cu02-planes of the high Tc cuprates. 

The phase diagram with t' /t = —0.7 and Ui/{At) — —0.5, and f/o/(4i) = is displayed in Fig. |5[ Due to breaking 
of electron-hole symmetry, the diagram is now a-symmetric around half occupation of the band. Roughly speaking 
s*-pairing is favoured far away from half filling of the band, whereas d-wave pairing becomes the most stable solution 
near half filling. We also notice from this plot, that the a-symmetry implies that the highest T^s and d-pairing 
superconductor are to be expected on the left-hand ('hole-doped') side of half-filling. Lower T^s and s-pairing occur 
on the right-hand side. 

Let us now consider the A/Tc-ratio following from the gap equation. Within the context of BCS theory we have 
Ao(r) = at Tc, so that Tc follows from 

=^eq:^tanh(^^^)[cos9^a±cosqj,a]^ (32) 

where the ± sign refers again to the two symmetries of pairing. This equation can be easily solved numerically. The 
result is, that for extended s-wave pairing the ratio 2Ao/fcsTc is 6.5, whereas for d-wave pairing it rises gradually 
from 4 if \Ui\ <C W, up to 6.5 in the limit where \Ui\ ^ W. This is not sensitive to the value of the parameter t' . We 
should keep in mind here, that Aq is the maximum value reached by A(fc) (respectively at the (tt, 0)- and (tt, 7r)-point 
for d- and s*-pairing). 

Let us now look how the mean field estimate of Tc depends on the coupling strength In Fig. || Tc'^^ /W is 

displayed as a function of |C/i|/VF for the d-wave channel. First of all we notice, that for |C/i| > W/A the value of T*^^ 
is about |J7i|/4. For << 1 this crosses over to a quadratic dependency T^^ = A\Ui\^ jW . For comparison 

a similar curve is displayed for conventional s-wave pairing, using the negative U Hubbard model in a band with a 
square DOS. We notice that the mean field transition temperature with the latter model becomes T^^ — |?7o|/4 for 
large \Uq\ (which is actually outside the range of validity of the BCS weak coupling approach |^,^), and has the 
familiar BCS-like exp (— VK/|L/o|) behaviour for small Uq. The Tc for the extended s-wave pairing lies again below the 
negative Uq curve, and is only finite above a threshold value of |J7i| as discussed above. 

Let us finally turn to the collective mode spectrum. We can anticipate, that again d-wave phase fiuctuations exist below 
the particle-hole continuum. In addition, because there is an on-site repulsive Uq, a branch of spin-fluctuations can be 
pulled below the particle-hole continuum. In Fig. ^the collective mode spectrum is displayed, using [/i/(4t) = —0.5, 
and rig = 0.85, and with C/o/(4i) ranging from to 1.5. In the plot for [/q = we already notice, that the particle-hole 
continuum has 8 points in /c-space where it touches the horizontal axis: The Fermi surface crosses the node-lines 
kx = ifcy at the coordinates (±(7r — 5)/2,±(7r — ^)/2), hence the particle hole spectrum is gapless for the Q- vectors 
(0, ±(7r — 5)), (±(7r — (5), 0) and (±(7r — 5), ±(7r — 5)). Precisely for these Q-values the spin (and charge) susceptibility 
acquires the largest value, also in the superconducting state, hence if we switch on a finite value of the repulsive 
on-site C/q, a spin-density wave starts to develope around the (±7r, ±7r) points on the Fermi surface. Clearly the 
ground-state is no longer of the form of Eq. |^, and the corrections may become strong enough to completely destroy 
superconductivity. As, on the other hand, the spin-density wave exists around a portion of the Fermi surface where 
the gap is zero (and therefore contributes the least to the ground-state energy), whereas the maximum gap- value is 
at the [±7r, 0] and [0, ±7r] points, there may actually be a coexistance of superconductivity and a spin-density wave 
in different portions of the Fermi surface. 

From Fig. |^ we can see, that the region taking part in the formation of the spin-density wave quickly spreads around 
the (±(7r — (5)/2,±(7r — 5)/2) points if C/o/(4t) increases, leaving a small region around [±7r,0] and [0, ±7r] for the 
formation of a superconducting condensate if ?7o/(4<) = 1. The phasediagram for [/o/(4t) = 1, and Ui/{At) = —0.5 is 
indicated in Fig. ^p. The shaded area roughly indicates the region with an instability towards a SDW. In principle 
a mixed SDW-superconducting state may exist for all concentrations. It is not possible to decide from the numerical 
results presented above whether or not there is a sharp phase boundary separating regions with a magnetic instability 
from superconducting regions. 
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V. CONCLUSIONS 



A unified approach is presented to the calculation of the collective modes of spin, charge, phase and amplitude 

in superconductors with a non-trivial pairing interaction. The expressions for the dynamical spin- and charge- 
susceptibilities are generalized to take into account superconductivity at general values of momentum and frequency. 
Several examples are treated. Notably the response functions of a layered charged electron gas, with a pairing inter- 
action in the rf-wavc channel is considered in the absence and presence of an on-site Hubbard repulsive interaction. 
An incipient instability toward a spin-density wave follows from the softening of the collective mode spectrum near 
Q = (tt, tt) in the d-wave paired state. 
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FIG. 1. Diagrams taken into account in the RPA. Exchange self energy (a), particle-hole mixing (a'), polarization vertex (b 
and b'), exchange scattering (c and c'), and direct particle-particle scattering (d and d'). Diagrams b', c' and d' exist only in 
the superconducting state. 

FIG. 2. Collective mode spectrum of a superconducting layered electron gas, assuming s-wave pairing. The parameters 
are: Ep/{4:t) = 0.35, Ry /{4:t) — 4.0. Q±c is varied with to tt with increments of 0.27r (top to bottom solid curves). The 
dashed curves are the boundaries of the region of Landau damping, (a) normal metal, U —0 and (b) superconducting state, 
U/ (At) = —0.67. (c) The amount of p (solid) and (p (dashed) character of the collective modes in Fig. |^ as a function of 
collective mode energy. The interruption occur where the modes become Landau damped. 



FIG. 3. Phase diagram in the Ui-n plane, where n is the number of electrons per unit cell, with t' = —0.7 and t/o/(4f) = 0. 

FIG. 4. Phase-fluctuating collective mode versus momentum for a layered electron gas with long range Coulomb interactions 
{Ry* — 20), an on-site repulsive interaction (Uo/W = 0.5), and a nearest neighbour attractive interaction Ui/W = —0.5. The 
number of electrons is We = 0.2 (a), rie = 0.25 (b), rie = 0.3 (c) and rie = 0.4 (d). 



FIG. 5. (a) Phase diagram in the T-n plane, where n is the number of electrons per unit cell, with t' = —0.7t and Uo / (it) — 0, 
and ;7i/(4t) = -0.5. (b) The same with f/o/(4t) = 1 

FIG. 6. Solid curve: Tc/\Ui \ calculated for the d-wave channel of the exchange-only model with t' = and 1 electron per 
site. The same curve is obtained for t' — —0.7t with 0.7 electron per site. Open lozenges: Tc of the s*-wave channel with the 
latter parameters. Dotted curve: Tc/\U\ versus |(/|/iy for the negative U Hubbard model taking a square DOS. 



FIG. 7. The collective modes in the d-wave paired state, using Ui/i/U) = —0.5, and — 0.85, and with C/o/(4t) =0 (a), 0.5 
(b), 1 (c), and 1.5 (d). 



14 



